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Analytical expressions are derived for sums of matrix elements and their squared moduli over 
many-body states with given total spin — the states built from spin and spatial wavefunctions 
belonging to multidimensional irreducible representations of the symmetric group, unless the total 
spin has the maximal allowed value. For spin-dependent one-body interactions with external fields 
and spin-independent two-body ones between the particles, the sum dependence on the many-body 
states is given by universal factors, which are independent of the interaction details and Hamiltonians 
of non-interacting particles. The sum rules are applied to perturbative analysis of energy spectra. 
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INTRODUCTION 

Calculations of quantum-mechanical system properties 
require matrix elements between its states. For complex 
systems, even a calculation of the matrix elements can 
constitute a complicated problem. However, in certain 
cases, sum rules can be derived from general principles, 
providing analytical expressions for sums of matrix ele¬ 
ments or their products. The Thomas-Reiche-Kuhn and 
the Bethe sum rules were obtained at early years of quan¬ 
tum mechanics. These and similar rules (see [1]) are for¬ 
mulated for weighted sums of oscillator strengths, which 
are proportional to squared moduli of transition matrix 
elements, over certain sets of eigenstates. The rules were 
applied to radiative transitions and scattering problems. 
Sum rules for dynamic structure factors (see [2]) are em¬ 
ployed to obtain information on collective behavior of 
many-body systems. Various sum rules are also used in 
nuclear and solid-state physics, as well as in quantum 
field theory. 

The present work derives sum rules for many-body sys¬ 
tems of indistinguishable spinor particles. The particles 
can be composite, e.g., atoms or molecules, and the spin 
can be either a real angular momentum of the particle 
or a formal spin, whose projections are attributed to the 
particle’s internal states (e.g. hyperfine states of atoms). 
In the latter case, the particle spin i means that only 
two internal states are present in the system. This for¬ 
mal spin is not related to the real, physical, spin of the 
particles, which can be either bosons or fermions. 

Many-body states of spinor particles can be described 
in two ways (see [3]). In the first one, each particle is 
characterized by its spin projection and coordinate, and 
the total wavefunction is symmetrized for bosons or an¬ 
tisymmetrized for fermions over permutations of all par¬ 
ticles [see Eq. (19) in Sec. I below]. The second ap¬ 
proach is based on collective spin and spatial wavefunc¬ 
tions. These wavefunctions depend on spins or coordi¬ 
nates, respectively, of all particles and form representa¬ 
tions of the symmetric group (see [4-7]). In the case of 
spin-i particles, the representation is unambiguously re¬ 


lated to the total many-body spin. If the total spin is less 
than the maximal allowed one {N/2 for N particles), the 
wavefunctions belong to multidimensional, non-Abelian, 
irreducible representations of the symmetric group (see 
[4-7]), beyond the conventional paradigm of symmetric- 
antisymmetric functions. The symmetric or antisymmet¬ 
ric total wavefunction — such functions only are allowed 
by the Pauli exclusion principle — is represented as a 
sum of products of the spin and spatial wavefunctions 
[see Eq. (2) in Sec. I below]. 

In the case of non-interacting particles in spin- 
independent potentials, all states with the given set of 
spatial quantum numbers are energy-degenerate and the 
two kinds of wavefunctions are applicable, related by a 
linear transformation. The effect of spin-independent in¬ 
teractions between particles was analyzed by Heitler [8] 
using the theory of the symmetric group irreducible rep¬ 
resentations. That work demonstrates that the average 
energy of states within given irreducible representation is 
proportional to a certain sum of the representation char¬ 
acters. The character dependence on the representation 
lifts the degeneracy of states related to different repre¬ 
sentations, and the wavefunctions with defined individ¬ 
ual spin projections become inapplicable. It is a gener¬ 
alization of the well-known energy splitting between the 
singlet and triplet states in two-electron problems. 

Although the derivation [8], being done at early years 
of quantum mechanics, did not take into account spin 
degrees of freedom and supposed that total wavefunc¬ 
tions can have arbitrary permutation symmetry, the re¬ 
sults remain valid for symmetric or antisymmetric to¬ 
tal wavefunctions, composed from spin and spatial func¬ 
tions of arbitrary symmetry. Matrix elements of spin- 
independent Hamiltonians between the latter wavefunc¬ 
tions can be reduced to the matrix elements between 
spatial wavefunctions due to orthogonality of the spin 
wavefunctions (see Sec. HI). Besides justification of the 
Heitler results, this reduction provides basis for spin-free 
quantum chemistry (see [6, 7]) — the method of cal¬ 
culations of energies and other properties of atoms and 
molecules. 

Spinor quantum gases are intensively studied start- 
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ing from the first experimental [9, 10] and theoretical 
[11, 12] works (see book [2], reviews [13, 14] and ref¬ 
erences therein). The collective spin and spatial wave- 
functions were used in derivation of exact quantum solu¬ 
tions for one-dimensional homogeneous gas [15, 16] and in 
analyses of selection rules and correlations [17]. SU{M)- 
symmetric gases, introduced in Refs. [18-20] and recently 
observed in Refs. [21, 22], are described in similar way 
[19], where the total wavefunction is composed of spin 
and electronic functions. 


Other forms of many-body wavefunctions with dehned 
total spin have been employed as well. The Lieb-Mattis 
theorem for ordering of energy levels in fermionic systems 
has been derived in Ref. [23]. One-dimensional gas of 
spin-i fermions in arbitrary potential has been analyzed 
for hardcore zero-range interactions in Ref. [24], where 
an exact solution was derived, and for zero-range interac¬ 
tions of arbitrary strength in Ref. [25], where qualitative 
properties of energy spectra are presented. An exact so¬ 
lution for one-dimensional hardcore Bose-Fermi mixture 
was derived in Ref. [26]. Intersystem degeneracies in 
spin-i Fermi gases and energy spectra for certain few- 
body systems have been obtained in Ref. [27]. Symme¬ 
tries of trapped and interacting bosons and fermions and 
qualitative behavior of the energy spectra at intermediate 
interaction strengths were analyzed in Refs. [28, 29]. 


The sum of matrix elements of spin-independent in¬ 
terparticle interactions directly follows from the Heitler 
results [8]. The present paper provides the sums of 
squared moduli of these matrix elements, as well as sums 
of matrix elements and their squared moduli for spin- 
dependent external fields. Such fields can be used for 
transfer of population between states with different total 
spins, as described in [17]. Besides, spin-changing ma¬ 
trix elements can provide an estimate of stability of the 
well-defined-spin states. 


Section I sets the analyzed problem and provides rep¬ 
resentations of spin, spatial, and total wavefunctions for 
separable spin and spatial degrees of freedom and for non¬ 
interacting particles. Wavefunctions with defined par¬ 
ticle spin projections are discussed in this section too. 
Section II contains derivation of the sum rules. Matrix 
elements of spin-dependent external fields for different to¬ 
tal spin projections are related using the Wigner-Eckart 
theorem. Then sums of these matrix elements and their 
squared moduli are calculated for the maximal allowed 
spin projections. Sum rules for spin-independent interac¬ 
tions between particles are provided in Sec. III. The sum 
rules are applied to description of the shifts and splittings 
of energy levels in Sec. IV. The quantitative properties 
of energy spectra are provided for arbitrary number of 
particles in the regime of weak interactions using pertur¬ 
bation theory. Appendix contains calculation of sums, 
used in Sec. II. 


I. THE HAMILTONIAN AND 
WAVEFUNCTIONS 

Consider a system of N particles with the Hamiltonian 

H = Hspat + 77spin (1) 

being a sum of the spin-independent Hspat and 
coordinate-independent Tfspin- Each of i7spat and Tfgpin 
is permutation-invariant. 

The total wavefunction is expressed in the form 

( 2 ) 

t 

rg\ ^ _ ,(S) 

Here spatial and spin ' functions form bases of 
irreducible representations of the symmetric group Sn of 
V-symbol permutations[4-7]. This means that a permn- 
tation V of the particles transforms each function to a 
linear combination of functions in the same representa¬ 
tion, 

P$if)=sgn(P)^Z4i^l(P)cI>(f2 

t' 

Here the factor sgn(P) is the permutation parity for 
fermions and sgn(P) = 1 for bosons. This factor provides 
the proper permutation symmetry of the total wavefunc¬ 
tion 


=sgn(P)<f (3) 

The matrices of the Young orthogonal representation 
[4-7] Dp\{V) of the symmetric group Sn are associated 
with the two-row Young diagrams A = [N/2+S,N/2 — S], 
which are unambiguously related to the total spin S. Dif¬ 
ferent representations, associated with the same Young 
diagram, are labeled by multi-indices n and I for the 
spatial and spin functions, respectively. The represen¬ 
tation basis functions are labeled by the standard Young 
tableaux t and t' of the shape A. The dimension of 
the representation is equal to the number of different 
tableaux. 


V!(2^+I) 

“ {N/2 -k A -k l)!(V/2 - S')!' 


( 4 ) 


If S = N/2, fs = I, Dpl(/P) = 1, and the functions 

„(S) 

and ^ remain unchanged on permutations of parti- 
cles or change their sign for fermions). Otherwise, 

the functions belong to multidimensional, non-Abelian 
irreducible representations of the symmetric group. For 
example, the states of IV = 3 particles with S = 1/2 are 
associated with the Young diagram [2,1] and there are 
fs = 2 standard Young tableaux with the Yamanouchi 
symbols (see [4, 5]) (2,1,1) and (1, 2,1). 
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The Young orthogonal matrices obey the orthogonality 
relation[6, 7] 

^ D^^'Xv)D^^X'P) = ^Su-Srr'Sxx,, ( 5 ) 

the general relation for representation matrices 

^ = D^^liVQ), (6) 

t 

and the relation for orthogonal matrices 

(7) 

Additional relations can be obtained for elements of 
the first column Dl^Qji'P) of the Young orthogonal ma¬ 
trices. Here [0] is the first Young tableau, in which the 
symbols are arranged by rows in the sequence of natural 
numbers. For example, the Young tableaux [0] have Ya- 
manouchi symbols (2,1,1), (2, 2,1,1), and (2,1,1,1) for 
the Young diagrams [2,1], [2^], and [3,1], respectively. 
Each permutation involving symbols between jmin and 
jmax can be written as a product of elementary trans¬ 
positions Vjj+I with jmin < j < jmax (see [5-7]). Ac¬ 
cording to the Young orthogonal matrix calculation rules 
(see [5-7]), Dl^{Vjj+i) = 6rt if j and j -I- 1 are in the 
same row of the Young tableau t. Then Eq. (6) leads to 
Zl|j[j](P) = jt[o] if file permutation V involves the sym¬ 
bols in one row only and can be, therefore, written as a 
product of elementary transpositions of symbols in the 
same row. Let V and V" be, respectively, arbitrary per¬ 
mutations of the symbols in the first and in the second 
row of the Young tableau [0], which do not permute sym¬ 
bols between the rows. Then we get, using Eq. (6), 

r 

( 8 ) 

The spatial and spin wavefunctions form orthonormal 
basis sets, 

(9) 

= Ss'sSt'tSin ( 10 ) 

The spatial functions of non-interacting particles are 
expressed as [6, 7] 

/ r \ 1/2 N 

\ -y -p j^i 

ai) 

in terms of the spatial orbitals — the eigenfunctions 
(/?ri(r) of the one-body Hamiltonian of non-interacting 
particle HqU), 

HoiA^Pnirj) = e„(/?„(rj), (12) 

where rj is the H-dimensional coordinate of jth parti¬ 
cle {D can be either 1, 2, or 3 in real physical systems). 


The representation is determined by the set of the spatial 
quantum numbers {n} and the Young tableau r, which 
can take one of fs values. Then the multi-index n is 
specifically chosen as r{n}. All quantum numbers rij in 
the set {n} are supposed to be different. This situation 
takes place in non-degenerate gases, when probabilities of 
multiple occupation of spatial states are negligibly small, 
although the multiple occupation is not forbidden by it¬ 
self. Another example is an optical lattice in the unit¬ 
filling regime [17]. 

The functions (11) satisfy the Schrodinger equation 

N N 

i=i i=i 

Their eigenenergies are independent of r. Therefore, 
there are fs degenerate states of non-interacting parti¬ 
cles for each set {n}. Tilde denotes wavefunctions cor¬ 
responding to the spatial Hamiltonian without interac¬ 
tions between particles. Then Eq. (2) gives us the total 
wavefunctions of particles with no coordinate-dependent 
interactions 

t 

In the absence of interactions between spins, the spin 
wavefunction are eigenfunctions of the total spin projec¬ 
tion operator Sz and can be expressed as 

N/2+S^ N 

■P i=l j=N/2+S^ + l 

(14) 

Here the multi-index I is specifically chosen as the total 
spin projection Sz- In the case of the spin wavefunction, 
each of two spin states, | f) and | j,), has to be occupied 
by several particles, if A^ > 2. Hence the normalization 
factor [30] 

^ ^_ 1 _ / i2S+l)iS + Sz)l 

{N/2 + Sz)KN/2 - A)! Y {N/2 + S+ 1)(25)!(^ - ^.)! 

(15) 

differs from the one in the spatial wavefunction (11). Be¬ 
sides, the Young tableau r can take now only the value 
of [0]. As a result, only one representation is associated 
with given total spin S and its projection Sz- The total 
wavefunction with the defined Sz is then expressed as 

( 16 ) 

t 

In combination with the spatial wavefunction (11), the 
spin wavefunctions lead to the total wavefunctions of 
non-interacting particles, 
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(again, tilde denotes that the wavefunctions involve spa¬ 
tial orbitals 93„(r) of non-interacting particles). There 
are fs wavefunctions, labeled by the Young tableau r, 
having the total spin S and the set of spatial quantum 
numbers {n}. Then the total number of wavefunctions 
with the given total spin projection Sz will be 

N/2 

fs = N\/[iN/2 + Sz)liN/2 - 5,)!]. (18) 

s=s^ 

In the alternative approach, mentioned in Introduc¬ 
tion, each particle has a given spin projection and the 
total many-body wavefunction is represented as (see [3]) 

N 

= {Nl)~^/^'^sgn{V)Y[(pn^{r-pj)\crj{Vj)), 

V j=i 

(19) 

where the spin projection aj can be either f or j, and 
given total spin projection Sz, the set {a} contains 
N/2-\- Sz spins f and N/2 — Sz spins j,. For a fixed set of 
spatial quantum numbers {n}, the number of such states 
is the number of distinct choices of N/2 + Sz particles 

with spin up, and is then equal to the number (18) of the 

~ / 0'\ 

states (17). Then the sets of degenerate states 

and ^'{n}{CT} can be related by an unitarily transforma¬ 
tion. For interacting particles, the energy degeneracy of 
states is lifted, as shown by Heitler [8] and will be 

discussed in Sec. IV, and such transformation becomes 
impossible. 


II. SUM RULES FOR ONE-BODY 
INTERACTIONS 


TABLE I. Coefficients in Eq. (23) 


k 

s- 

V 


0 

1 

-1 

V(S-S^-I-1)(S-I-S*) 

/(S-I-S^-1)(S+S^) 

V2S 

V 25(25-1) 

0 

Sz 

/ S2-S2 

s 

V S(2S-1) 

1 


/(S-S;,-1)(S-S*) 

V2S 

V 25(25-1) 


are the spin raising and lowering operators for jth par¬ 
ticle. The interaction Uq conserves the z-projection of 
the total many-body spin, while U±i raises or lowers it. 
The interaction of the spin-up or spin-down state can be 
expressed in terms of Ug and the scalar U, 

17t = ^t/(r,)|t(j))(t {j)\=Uo + lu 

' ( 22 ) 

U^ = Y^ U{r,)\ i (j))U (J)l = -Uo + ^u. 

3 

Consider matrix elements of the spherical vector and 
scalar interactions between eigenfunctions (16) of Sz- 
Their dependence on Sz follows from the Wigner-Eckart 
theorem (see [32]). The matrix elements of the spherical 
scalar (20) are diagonal in spins and independent of the 
spin projection, 

According to the Wigner-Eckart theorem, the matrix el¬ 
ements of the spherical vector components (21) can be 
factorized into the 3j-Wigner symbols and the reduced 
matrix elements 


A. The spin-projection dependence 




iS) 


\S'-S' 


S' 1 S 
-S'' k Sz 


(n',S'||U||n,S). 


Permutation-invariant interactions of particles with ex¬ 
ternal fields can be expressed in terms of the spherical 
scalar 


U = Y^U(v,) (20) 

3 

and three spherical vector components 

Co = ^ U{v,)sz{3), U±i = E U{rj)s±{j) 

( 21 ) 

(see [31]). Here 

= ^(11 (j))(t (j)|-U(j))U (j)l) 


Then the reduced matrix elements are expressed in terms 
of the matrix elements oiUk for the maximal allowed spin 
projection 

(n',S'||C||n,S) = ( ^ ^ 


and the matrix elements with arbitrary spin projections 
can be expressed as 


,(A) lU. i,T,(S) \ _ 


{Ky'\Uk\Kl)=Ss'^s.+kX^s. 


nS ) 

(23) 


with the factors 


X 


{S,S',q) 

S,,k 




S S' q\ 

Sz -Sz-k k ) 


(s S' q \ 
\S -S' S' - S J 


is the z-component of the spin and 

5+(j) = it(j))a(j)i, s_(j) = u(j))(t(j)i 


X 


-1 
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Here S' < S and, according to the properties of the 
3j-Wigner symbols, the matrix elements (23) vanish if 
[S' — S"| > 1 (in agreement to the selection rules [17]). 
Values of non-vanishing coefficients, calculated with the 
3j-Wigner symbols [3, 32], are presented in Tab. I. Her- 
mitian conjugate of Eq. (23), together with relations 
U+i = —U-i and Uq = Uq, gives us the matrix elements 
for S' = S+1. 

Thus, each matrix element of a spin-dependent one- 
body interaction with an external field is related to ma¬ 
trix elements for the maximal allowed spin projections, 
which will be evaluated in the next section. 


B. Matrix elements for non-interacting particles 


Matrix elements of the spherical scalar (20) can be 
evaluated exactly for general spin wavefunctions. Due to 
the orthogonality of the spin wavefunctions (10), the ma¬ 
trix elements are diagonal in spin quantum numbers and 
can be reduced to the matrix elements between spatial 
wavefunctions, 

(24) 

t i 

Let us calculate the spatial matrix element for the general 
case, S S', having in mind further analysis of spherical 
vectors. Equations (11) and (7) lead to 


\ t'r'li 




Vfsfs' 

m 


^sgn(Q)DM(Q) 

'R.,a 


X sgn(7^)i:>tt^(7^)((/?„/gJC/(r,)|v5„^,) 

i' 


the identity permutation, one finally gets 

N 

= 5ss'5w5r'rY.^n'^\U\n,) H 

i=i 

(26) 

It is a special case of the matrix elements obtained by 
Heitler [8] and Kaplan [6]. 

For the spherical vector interactions (21), the matrix 
elements cannot be represented in so simple a form. How¬ 
ever, rather simple expressions can be derived for sums 
and sums of squared moduli of the matrix elements be¬ 
tween eigenfunctions of Sz- It is enough to consider ma¬ 
trix elements of [/_i and the spin-up state interaction 
[Tf- for the maximal allowed spin projection, S'' = S', 
Sz = S, as Eq. (22) and the Wigner-Eckart theorem (23) 
relate to them each matrix element of each interaction. 
In the basis of the non-interacting particle wavefunctions 
(17), the matrix elements of l/f can be decomposed into 
the spatial and spin parts, 

X t (*))(t (*)|Si^^ 

(27) 

The spatial matrix elements are given by Eq. (25). The 
spin matrix elements include projections of the spin wave- 
functions (14) 


(t (*)|s|^))=C'ssE<d](^)E'^*^' 

V 1^1 

Ai N 

xnit(^7)) n u(^4))- 

j=M+i 


Substituting V = QVix-iWe get 


The Kronecker 5-symbols appear here due to the orthog¬ 
onality of the spatial orbitals (pn and the absence of equal 
quantum numbers in each of the sets {n} and {n'}. Due 
to the 5-symbols, all but one spatial quantum numbers 
remain unchanging. Supposing that the unchanged nj' 
are in the same positions in the sets {n} and {n'}, one 
can see that the Kronecker symbols lead to Q = 77., and, 
therefore, 




Vfsfs' 

m 


TZ 


'X {n'm\U\nnz) H 

j'yi'R.i 


where {n'\U\n) = f d^rtp^, (r)U(r)(p„(r). Then, substi¬ 
tuting this expression into (24), using (6), (7), and the 
property of representations = Sr'r, where £ is 


(t (*)|S||)) = CssEE^1[c1](2^'>i)'^*cx, 

Q /-I 

Ai-1 N 

Xnit(27)) n It(27)). 

7=1 7=>'i-l-i 

The permutation Vix^ permute symbols in the first row 
of the Young tableau [0]. Therefore, D||[jj(Q7^;Ai) = 

D|j[|j(Q) [see Eq. (8)] , the summand in the equation 
above is independent of I, and the projection can be ex¬ 
pressed as 

(t (*)|Si^^)=AiCssE<cl](2)<5iQA, 

Q 

Ai-1 N 

X nit (27)) n lt(27)). 

7=1 7=Ai+i 
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The projection involved into matrix elements of U-i is 
evaluated in the same way, 

(i = A2Cs5E^lfd](2)'^*C(A, + l) 

s 

Ai N 

xnit(2j)) n \iiQj))- 

j = l j=Ai+2 


where A' = [Ai — 1, A 2 + 1], is calculated in a similar way. 

The explicit expressions (28) and (29) are rather com¬ 
plicated as they include Young orthogonal matrices and 
summation over all permutations. The next section pro¬ 
vides expressions for sums of the matrix elements and 
their squared moduli, which are much simpler. 

C. Sum rules 


In the spin matrix elements of U-f, 

(4s'l t (*))(t (*)|sif^) = Sss' [AiC'ss]"^^lJl](2)<5*eA, 

Q 

Sn,QV'V" 

n 'P','P" 

the permutations TZ and Q can be different by permuta¬ 
tions V' of the first Ai — 1 symbols and V" of the last 
A 2 ones. As the permutations V' and V" do not per¬ 
mute symbols between rows in the Young tableau [0], we 
have d|^jpj( 7^) = D|^jgj(Q) [see Eq. (8)]. Since the num¬ 
bers of permutations V and V" are (Ai — 1)! and A 2 !, 
respectively, the spin matrix elements take the form, 

(S|fil t (*))(t (*)|si^^) = (Ai - 1)!A2 !A?C'|s 

Q 

Let us substitute this equation and (25) into (27), per¬ 
form the summation over t and t' using Eq. (6), and 
substitute V = j = TZi. Then the Kronecker 

symbol leads to Vj = Q~^i = Ai, and we get 

V 

N 

X ^(5AiT’j(n'|C/|nj) (28) 

i=i jYj 

The matrix element 


The sum of diagonal in total spin S and r matrix ele¬ 
ments can be written out as 


r 

N 

X n 

i=i i'/i 

The universal factors are independent of the 

matrix elements (n'|17|nj). For U-f, the factor 
Y('^)[{7^] can be derived from (28) using the equalities 
Er-D[ o]U^)^[o] 1(^) = -^[o][o](^) = 1 [obtained with (6) 
and (7) [ and = (n — 1)!, as 

Y(®)[C/t] = Y+'5'- (30b) 

It is equal to the number of the spin-up atoms. For the 
spherical vector component C/q, the factor Y(‘^)[[/o] is ob¬ 
tained using Eq. (22), 

Y^^'>[Uo] = S. (30c) 

Equation (26) leads to 


(^lE'Is-il^-il^Sls) = ^Ai!A2!(A2+1)CssCs-is-i 

N 

XE ^[oh'(^)^[0]l i'P) E in' \U\n,)l[ , 


y(^)[I7[ = N. 


(30d) 


i=i 


t'/i 


The sum of squared moduli of the matrix elements (28) 
(^®) and (29) can be expressed, using Eqs. (4) and (15), as 





0"¥^3 


n E 




(31) 
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(32) 




where 




-)[^'] /'■E>^n[^l 


[-^'1 (n\ nt-^l 


(33) 


r,r' V 


Q 


These sums are calculated in Appendix. It is shown that 


(S,S) _ mjN-iy. r _ Az 

" fsXl [ ' Ai - A 2 + 2 

^(s-i,s) _ N\{N - 1)! 

'' “ /sAi 

are independent of j , and 


(34) 

(35) 


y,{S' fS) 


N\{N -2)\ 
fs 


Sss' - 


1 


yy 

N-1 


iS',S) 


(36) 


for any j' ^ j. 

If the sets of spatial quantum numbers {n} and {n'} 
are different, the product of Kronecker symbols in (31) 
and (32) does not vanish only if j = j'. Then the sum 
of squared moduli of the matrix elements can be written 
out as 


r r' 


|f>a|4'. 




N 


N 


i=i i'/j 

where S' < S and the difference S — S' is unambigu¬ 
ously determined by the operator Ua- Each term in the 
sum here changes one spatial quantum number, conserv¬ 
ing other ones. If U{y) = const, the sums vanish since 
{ypn'\U\ipn) = U{(pn'\ypn) = 0 for Ti ^ n'. The universal 
factors Y^^’^yUa, Ua], which are independent of the ma¬ 
trix elements {n'j\U\nj), are expressed in terms of Sjj 
Then Eqs. (34) and (35) lead to 


N 


S- 


N-2S 


2 ' ■ 4(5'-bl) 

- 25 -b 2 


(37b) 

(37c) 


‘± 

and Eq. (26) gives 

Y‘^^’^\U,U] = N. 

The factor {/p] for the spherical vector compo¬ 

nent Uq is obtained using (22). Since the matrix elements 
of IJ are diagonal in r [see Eq. (26)], one gets 




(37d) 


For transitions conserving the spatial quantum num¬ 
bers, {n'} = {n} and the Kronecker symbols in (31) and 
(32) are equal to one for any j and j'. Then sums of 
squared moduli of the matrix elements can be represented 
as 




Y,^^'°\Ua,Ua]{Uf 


+ r/^’°^[17a,17j(AC/)2j (38a) 


where 


(U) = ^j2<nyU\n,) 

i=i 


is the average matrix element and 


(At/) 




is the average deviation of the matrix elements of U{r). 
The universal factors YQ^’^\Ua,Ua] and Y}^’^\Ua,Ua] 
are independent of the matrix elements {nj\U\nj). Equa¬ 
tions (34), (35), and (36) lead to 

Yo^^’°\Ua, Ua] = (38b) 


(where defined y*-‘^^[t/_i] = 0, in addition to Eq. (30)), 
and 


yf’°^[i7_i,t7_i] = 


A(A- 25-b2) 


Y 


4(N-1) 

(^’°^[t7t,t7t] = r/^’°^[t7o,17o] 

S(N-2S)(N + 2S + 2) 
4(S + 1)(N-1) ■ 


(38c) 


(38d) 


If U(r) = const. At/ = 0, and, therefore, 

Er.r'= 0- Indeed, in this 
case, the spatial matrix elements (25) are equal to zero 
due to the orthogonality of the spatial wavefunctions with 
different spins. 

Thus, sums of matrix elements and their squared mod¬ 
uli are expressed in terms of universal factors, which are 
independent of the spatial orbitals and details of the ex¬ 
ternal fields, and sums of one-body matrix elements (or 



their squared moduli), which are independent of many- 
body spins. The sum rules, combined with the spin- 
projection dependence (23), provide information on each 
matrix element for an one-body spin-dependent interac¬ 
tion with an external field. 


III. SUM RULES FOR TWO-BODY 
SPIN-INDEPENDENT INTERACTIONS 

The permutation-invariant interaction between parti¬ 
cles is given by 

U=^U(r,-r,0. (39) 

i/i' 

Without loss of generality, we can restrict consideration 
to even potential functions, V{y) = U(— r), since their 
odd parts are canceled. Matrix elements of this inter¬ 
action can be evaluated for general spin wavefunctions. 
Due to the orthogonality of the spin wavefunctions (10), 
the matrix elements are diagonal in spin quantum num¬ 
bers and can be reduced to the matrix elements between 
spatial wavefunctions. 




$ 


( 5 ) V 
tr{n}l 


t i<.i' 


(40) 


(this reduction is used in spin-free quantum chemistry [6, 7]). Then, using (11), (39), and the property (7) of the 
Young orthogonal matrices, the spatial matrix elements can be expressed as 


= M E sgn(Q)DW(Q)sgn(7^)DW(7^) 


n,Q 




(41) 


i' z^i" 


The Kronecker 5-symbols appear here due to the orthogonality of the spatial orbitals (pn and the absence of equal 
quantum numbers in each of the sets {n} and {n'}. Due to the 5-symbols, all but two spatial quantum numbers 
remain unchanging. Supposing that the unchanged rii" are in the same positions in the sets {n} and {n'}, one can 
see that the Kronecker symbols allow only Q — TZ or Q = TZVa'. Then substitution of (41) into (40), using (6) and 
(7), leads to 




iS) 


TV! EE n 

TZ i<i' TZi'^j”^IZi 

5r'r{n'Ti,n'-^,,\V\nmnTZi') + sgn(T’^*/)D|.■)|.(7^T’^*/7^■^)(n)^,,r^)J,|U|r^7?,^r^7?,^')j , (42) 


X 


where {n[n 2 \V\nin 2 ) = f d^rid^r 2 p^,Jri)<p*^(r 2 )V(ri - r 2 )p„,(ri)<p„^(r 2 ). 
Taking into account that 


Wii'V ^ = VviVi' 


(43) 


(see [7]) and substituting TZi = j, one finally gets 




i<i' 


( 44 ) 


It is a special case of the matrix elements obtained by Heitler [8] and Kaplan [6]. 
The sum of diagonal elements of the representation matrix, the character 

xs{c) = J2Dl^Jir), 
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TABLE II. Characters Xs{C) of the classes C of conjugate elements of the symmetric group Sn of permutations of N symbols 
in the irreducible representations, corresponding to the spin S. The characters are calculated with the Frobenius formula[7, 34] 
and scaled to the representation dimension fs- 


c 

Xs{C)/fs 

{2} 

4S'’-l-JV'’+4S-4Af 

2iV(iV-l) 

{3} 

12S''+Af^ + 12S-10Ar 

4iV(iV-l) 

{4} 

N*-24N'-'+4N-‘(&S-‘+&S+29)-16N(10S'^+10S+9) + 16S(S+l)(S'‘‘ + S+12) 

8JV(]V-l)(JV-2)(]V-3) 

{2^ 

]V'‘-12]V'‘+8iV''(S'’-l-S-l-7)+8Af(10S'' + 10S-l-9) + 16S(S-l-l)(S''+S+6) 
4Ar(iV-l)(iV-2)(iV-3) 


is the same for all permutations V, which form the class of conjugate elements C [4-7]. Table II presents the characters 
for the classes appearing here. (Supplemental material for [17] contains a code based on the explicit expressions [33] 
for the characters.) The conjugated classes of the symmetric group Sn are characterized by the cyclic structure of the 
permutations. All permutations in the class C = {N'^^ ... 2'^^} have vi cycles of length /. This class notation omits 
\i vi = 0 and the number of cycles of the length one, i.e. the number of symbols which are not affected by the 
permutations in the class. This number is determined by the condition ~ ^■ Permutations of two symbols 

form the class {2}. This leads to the sum of diagonal in r matrix elements 




2 E ±Xs({2})(n',n'|y|njnj/)] Y[ 

i<i' 3'¥^3"¥^3 


(45a) 


where the sign + or — is taken for bosons or fermions, respectively. Similar expressions have been obtained for the 
total energy [8] and arbitrary observables [17]. If {n'} = {n}, the Kronecker symbols are equal to one for any j and 
f and the sum can be transformed to the form 

- l)/s ± • (45b) 


Here and above, the dependence on many-body states is given by universal functions fs and xs({2}), which are 
independent of the matrix elements (n]^n 2 |H|nin 2 ), while the average matrix elements 

2 _ 2 

^ N(N-l) ^ = N(N-l) ^ i'^3'^3\V\njnf) (46) 

of the direct and exchange interactions, respectively, are independent of the many-body states. 

Calculating the sum of squared moduli of the matrix elements (44), one can see that if the sets of spatial quantum 
numbers {n} and {n'} are different by two elements, the product of Kronecker symbols in the product of the matrix 
elements does not vanish only if the pair j, j' is the same in both matrix elements. Then the sum can be expressed as 


EK^''{n'pl^lnW/)l'=4/5E n 


I (n' n', I y I rij n/) I ^ -b I (n', n' | K | ) | ^ 

Xs({2}), 


± 2 - 


fs 


-Re ((n' n', Iy I Ujrif ) (n'/ n' |y | rij n^/)*) 


(47a) 


Here the equality Err' = Er-^ir {S) = fs was used. Each term in the sum above changes 
two of the spatial quantum numbers, conserving other ones. The case of a single changed quantum number will be 
considered elsewhere. 

For transitions conserving the spatial quantum numbers, {n'} = {n} and the Kronecker symbols in (44) are equal 
to one for any j and j'. Then 


r r*' 


fsN^N - i)^(y)E ± 2 xsm)N\N - i)^(y)dir(y)ex 

+ E EE 'b’lr' \y\nnnx^){nx^nj,\V\nj^nj^) 

iiAil h¥^32 '' 
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The trace of the Young matrix can be transformed in the following way (since ji ^ and j 2 ^ j' 2 ) 

)=xs({2n)+ ){xsm)-Xsi{2^})) 

r 

+ ^h32^j'ij2 ){fs-2xsm)+xs{m)), 

since G {3} for j[ ^ = S, and xs{£) = fs- Here and in what follows, Xs({3}) and xs({2^}) 

have to be equated to zero at iV < 3 and < 4, respectively, when the corresponding permutations do not exist. 
Using the identity 2/s + 4(A^ —2)xs({3}) + (iV — 2)(iV — 3)xs({2^}) = N{N— l)x|({2})//s (it can be directly proved 
with the characters in Table II), the sum of squared moduli of the matrix elements can be represented as 


with the universal factors 

y^(g.o) yj ^ _ ^^2 xs({3}) - xs({2 }) ^ = 2N{N -1) (l - ~ j . ( 47 c) 


fs 


fs 


Here 




(Ajr)^ = 


J=1 

2 


— E 

N -I ^ 


{nj>nj\V\njnf) - {V) 


—-— Y^{{nj,nj\V\n,n,,) - {V),^f 
^ ’ 3 < 3 ' 


(48) 

measure the average deviation of the exchange matrix 
elements. 

Thus, sums of matrix elements and their squared mod¬ 
uli are expressed in terms of universal factors, which are 
independent of the spatial orbitals and interaction po¬ 
tentials, and sums of two-body matrix elements (or their 
squared moduli), which are independent of the many- 
body spins. The universal factors are expressed in terms 
of characters of irreducible representations of the sym¬ 
metric group. The characters are functions of the total 
spin and the number of particles. 


IV. MULTIPLEX ENERGIES FOR 
WEAKLY-INTERACTING GASES 

As an example of applications of the sum rules, con¬ 
sider splitting of degenerate energy levels due to weak 
two-body spin-independent interactions. The Hamilto¬ 
nian of the system is a sum of one-body Hamiltonians 
HoU) of non-interacting particles and two-body interac¬ 
tions (39), 


N 


4pat =EE(j) + H 


(49) 


1=1 


The interactions split energies of the degenerate states 
(13). In the zero-order of the degenerate perturba- 


I- 

tion theory [3], the eigenenergies Esn (counted from the 
multiplet-independent energy of non-interacting particles 


X]j=i ^rij) are determined by the secular equation 

/ . ''rr' ^nr' ~ ^Sn^nr ) 


(50) 


(S) 

where Anr are the expansion coefficients of the wavefunc- 
tion (16) in terms of the wavefunctions of non-interacting 
particles (17), 


nSz / . rir ^ r{n 


n}Sz 


(51) 


and the matrix elements of the spin-independent two- 
body interaction (44) 




r{n}S- 


) 


do not couple states with different spins. 

Consider at first the case when the matrix elements 
Vdir = (uin 2 |U|nin 2 ) and 14x = (iiiii 2 |U|n 2 ni) are in¬ 
dependent of the spatial quantum numbers. E.g., this 
can take place in the case of zero-range interactions 
U(r) = U(l(r), if the spatial orbitals have a form of plane 
waves. In this case, the summation over TZ in the ma¬ 
trix element (42) for {n} = {n'} can be performed using 
Eqs. (6), (7), and the orthogonality relation (5) in the 
following way [5] 

E Dl^UnVwn-^) = E ) E (^) 


n 


TZ 


Nl 

— 'J^^r'rXsi{2}). 

Then the matrix elements become diagonal in r, 

Xs({2}) 


= 6rr'N{N - 1) (^Vdir ± 


fs 


- 14 , 
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FIG. 1. (Color online) Scaled average multiplet energies Es 
as functions of the multiplet spin S for N bosons or fermions. 
The energies are calculated with Eq. (52) and scaled charac¬ 
ters from Tab. II, assuming (I^}dir = (V) ex “ W- 



FIG. 2. (Color online) Scaled root-mean-square energy widths 
of multiplets as functions of the multiplet spin S for N parti¬ 
cles, calculated with Eq. (53). 


where the character y;5({2}) is given in Tab. II and the 
sign -|- or — is taken for bosons or fermions, respectively. 
The secular equation (50) is then satisfied by the eigen¬ 
vectors = Snr and eigenvalues Esn = VrP. Then 
all eigenstates with the given spin remain degenerate in 
energy. However, states with different total spins have 
different energies. 

In the general case, when the matrix elements of V 
depend on the spatial quantum numbers, the energies 
Esn can not be expressed in a simple form. However, 
using the equivalence of the sum of matrix eigenvalues 
to its trace and the sum of matrix elements (45b), the 
average multiplet energy can be expressed as 




1 

fs 




= N{N - 1) 



± 


X5({2}) 

fs 



(52) 


where the average interactions (H)dir and (4^)ex are de¬ 
fined by (46). (Here and below, the summation over n 
means the summation over states of interacting particles 
in a given spin-multiplet with a given set {n}.) It is 
a particular case of the general expression obtained by 
Heitler [8]. The average energies are plotted in Fig. 1. 

As the interaction lifts degeneracy of states with differ¬ 
ent total spins, transformation of the set of states with 
defined total spins to the set of states with given spin 
projections of particles becomes impossible. Then the 
former set remains the only valid set of eigenstates of 
interacting particles. 

For fermions, the average multiplet energy decreases 
with S. The Lieb-Mattis theorem [23] predicts opposite 
dependence. However, this theorem is formulated for the 
lowest-energy states with given 5, which can involve dif¬ 
ferent sets of {n} and have multiple occupation of spa¬ 


tial orbitals. In contrast, the average energies (52) are 
obtained for the fixed set of {n} and single occupations. 

The root-mean-square energy width of the spin-S' mul¬ 
tiplet (AEs) is defined by 

{AEsf ^ ^ - Esf = ^ E - ^5 


Due to orthogonality of the expansion coefficients, 
the secular equation (50) can be rewritten in the form 
Esn^n'n — leading tO 


1 

fs 


^ ^ \Esn^n'n\ 

n,n' 


1 

fs 


EKwKv. 


Then Eq. (47b) gives us 

{AEsf = TE°^[t>, H](AiF)" + 1>](A2'F)^ 

^ 53 ) 

where the universal factors [I/, V] and ^^] 

are expressed in terms of the representation characters 
by Eq. (47c), and the matrix element deviations (AiH) 
and {A 2 V) are defined by Eq. (48). The multiplet en¬ 
ergy widths are plotted in Fig. 2. If the matrix elements 
of V are independent of the spatial quantum numbers, 
(AiH) = {A 2 V) = 0 and, therefore, {AEs) = 0, in agree¬ 
ment with the above-mentioned degeneracy of states with 
given S in this case. The energy width is determined by 
characters, which were identified by Dirac [35] as con¬ 
stants of motions, corresponding to permutation symme¬ 
try, according to generalized Noether’s theorem. There¬ 
fore, the energy width can be considered as a conserved 
physical observable, related to this symmetry, as well as 
the average multiplet energy and correlations [17]. 

Using characters from Tab. H, the exact expression can 
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be approximated at » 1 by 
/V2 _ 45*2 

{^Esr « y^[2jV(4^^ - 3N){A^vy 

+ - AS‘^){A2Vf]. 

Consider now external fields described by one-body in¬ 
teractions. Matrix elements of a spin-independent field 
(26) are independent of r and spin quantum numbers. 
Therefore, this field leads to the same shift for all states, 
corresponding to the given set of spatial quantum num¬ 
bers {n}. In the first order of the perturbation the¬ 
ory, this shift will be ^f=i{nj\U\nj). Even strong spin- 
independent field leads to the same shift of all states, 
as it can be incorporated into the Hamiltonian of non¬ 
interacting particles. Then, the Schrodinger equation 
(12) will contain Ho{j) + U{rj). This leads to different 
one-body eigenfunctions (pn(j) and eigenvalues e„, but 
does not change the form of many-body wavefunctions. 

Spin-dependent spatially-homogeneous interactions 
[Eqs. (21) and (22) with U = const] commute with 
the spatial Hamiltonian of interacting particles (49). 
Since the spin wavefunctions (14) are eigenfunctions 

(S) 

of such interactions, the eigenfunctions 4')^^ of i?spat 
will be eigenfunctions of the Hamiltonian Hgpat + Uo- 
The energy shift of the states of non-interacting par¬ 
ticles due to the field Uq is equal to the matrix el¬ 
ement = Sr’rSzU. It is deter¬ 

mined by Eqs. (22), (23), (26) and (28), taking into 
account that {n'\U\n) = USnn’- The energy shift of the 
states of interacting particles (51) will be the same, as 
= S^UY^rKrAur = S.U. 

The spin-independent inhomogeneous and spin- 
dependent homogeneous fields, considered above, are 
consistent with the separation (1) of the spin and spatial 
Hamiltonians. If the external field depends both on spins 
and coordinates, this separation is violated, invalidating 
the use of collective spin and spatial wavefunctions for 
non-interacting particles. Nevertheless, these wavefunc¬ 
tions remain applicable to interacting particles whenever 
the external field is weak enough, and the energy shift 
can be estimated in the first order of the perturbation 
theory. The average shift is calculated using orthogo¬ 
nality of the coefficients Anr, Eqs. (23) and (30) in the 
following way, 

N N 

i=i i=i 

(54) 


CONCLUSIONS 

The symmetric group methods allow to evaluate the 
matrix elements of spin-dependent external fields (21) 
and spin-independent two-body interactions (39) in the 
basis with collective spin and spatial wavefunctions (2). 
These matrix elements agree to the selection rules [17]. 
For the matrix elements of spin-dependent external fields, 
explicit dependence on the total spin projection (23) is 
obtained using the Wigner-Eckart theorem. Analytical 
expressions are derived for sums of these matrix ele¬ 
ments (30) and their squared moduli [Eqs. (37) and (38)] 
over irreducible representations for both spin-conserving 
and spin-changing transitions. Dependence on the many- 
body states in these sums is given by the 3j Wigner sym¬ 
bols and the universal factors and 

These factors are independent of details of one- 
body Hamiltonians and external fields and are expressed 
in a rather simple form in terms of the total spin and 
number of particles. For spin-independent two-body in¬ 
teractions, the sums of matrix elements (45) and their 
squared moduli (47) depend on the many-body states 
only through the representation characters, which were 
identified by Dirac [35] as constants of motions, corre¬ 
sponding to permutation symmetry. The sum rules can 
be applied to the evaluation of energy-level shifts (54), 
splitting of states with different total spins (52), and 
spin-multiplet energy widths (53). Other possible ap¬ 
plications of the sum rules include estimates of the spin- 
multiplet depletion rates due to spin-dependent pertur¬ 
bations, as well as the population transfer rates between 
spin-multiplets using the spatially-homogeneous spin¬ 
changing and spatially-inhomogeneous spin-conserving 
pulses [17]. 


Appendix: Calculation of the sums (33) 

Using the relations (6) and (7) and substitution TZ — 
Q'P~^, the sum (33) can be represented in the following 
form 

n. V 

where A = [N/2 + S, N/2 - S] and A' = [iV/2 -f S', N/2 - 
S']. 

For j = j', there are {N — 1)! permutations V such 
that Vj = Ai. Then 

^ (A.l) 

n 

is independent of j. 
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For j ^ j', we have 


and I'l > Ai- Then [36] 


= X! '^Ai ,ni ^ Si^-pj'6\^^-pj 

V l^\i V 

= {N-2)\Y,K,ni{l-Ki) = {N-2)\{1-5xm). 

I 


Then 


^n' 


n 

_ N\{N -2)\ 

Is 


[A] 

[ 0 ][ 0 ] 

^AA' 


{n){N-2)\{l-5x,,nx,) 


_}_y(S',S) 

N - 1 ’ 


(A.2) 


where the last transformation uses Eqs. (5) and (A.l). 
The last expression in (A.2) is independent of j and j' 
and equivalent to (36). 



(- 1 )"' 


- nexV- 

Ai! 


Due to the Kronecker symbols in Eq. (A.l), the per¬ 
mutations V' do not affect Ai and zj, < Ai — 1. Therefore 
there are (Ai — 1)! permutations V' , A2! permutations 
P", and number of distinct choices of the sets of and z)! 
are given by the binomial coefficients ('''^~^) and (^^), 
respectively. Then for S = S' Eq. (A.l) can be trans¬ 
formed as follows, 


A, 


(S.S) 


= (1V-1)! (Ai-l)!A2!f^' 

r, \ ^ex J 

P^ex —0 


A2 

'^ex 


Ai 

‘^ex 


-2 


(A^ - 1)!(A2!)^ (Al-nea:)!,, A 


A? 


„.:lo (A2-ne.)! 


The sum over Uex can be calculated, leading to (34) . 
If S" = S' — 1 we have 


The Young orthogonal matrix elements in (A.l) have 
been calculated by Goddard [36] in the following way. 
Each permutation P can be represented as 




fc=l 


where P' are permutations of symbols in the first row of 
the Young tableau [0] (Ai Hrst symbols), P" are permu¬ 
tations of symbols in the second row (A 2 last symbols), 
and Pi'y'^ transpose symbols between the rows as zj, < Ai 


A, 


(S-l.S) 


jj 


= (»-!)! i;(A.-i)!w(^;-‘)(;^)( 

Uex^O ^ ^ V^ex/ \ 


'^ex 


Ai-l 

'^ex 


(iV-l)!(A2!)2 ^ (Ai-ne.)! 


E 


TL^x —0 


(A2 ZZea:)! 


giving (35) . 
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